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The electroluminescence of molecules confined inside a nanocavity in the scanning tunneling mi- 
croscopy possesses many intriguing but unexplained features. We present here a general theoretical 
approach based on the density matrix formalism to describe the electroluminescence from molecules 
near a metal surface induced by both electron tunneling and local surface plasmon excitations 
simultaneously. It reveals the underlying physical mechanism for the external bias dependent elec- 
troluminescence. The important role played by the local surface plasmon on the electroluminescence 
is highlighted. Calculations for porphyrin derivatives have reproduced corresponding experimental 
spectra and nicely explained the observed unusual large variation of emission spectral profiles. This 
general theoretical approach can find many applications in the design of molecular electronic and 
photonic devices. 



Light emission of molecules induced by the localized 
tunneling current in a scanning tunneling microscope 
(STM) has been studied for almost two decades |ll4l3l|. 
owing to its rich new physics and fascinating potential 
applications in molecular scale optics and electronics 



fplici 



14| . In principle, the electroluminescence of a 



molecule can be generated once an electron-hole pair is 
formed inside the molecule. This could happen when 
the highest occupied molecular orbital (HOMO) and the 
lowest unoccupied molecular orbital (LUMO) enter into 
the window of the external bias0, 01 • Such an electron 
driven mechanism has been explained by several groups 
using different theoretical methods]^, [ll|, [l^l. Very re- 
cently, a fascinating new experimental study by Dong et 
g/. [14| has shown that the electroluminescence spectra 
of tetraphcnyl porphyrin (TPP) molecules can be effec- 
tively modified by tuning the energy profile of the reso- 
nant nanocavity plasmons (NCPs) in a STM. The reso- 
nant hot electroluminescence arising directly from higher 
vibronic levels of the singlet excited state for porphyrin 
molecules, even upconversion luminescence, have been 
observed for the first time. It was suggested that the 
local nanocavity plasmons behave like a strong coherent 
optical source with a tunable energy However, there 
is no fully developed theoretical model ready to describe 
these new phcnomcnons. 

It was known that the light induced fluorescence and 
Raman scattering of molecules in the metallic nano- 
gap could be uniformly treated by the density matrix 
approach [l^ In this letter, we will show that the 
density matrix formulation can be generalized to describe 
the electroluminescence from molecules near a metal sur- 
face in a STM. In our model, the surface plasmons are 
regarded as strong coherent electromagnetic sources, the 
electron tunneling rates are determined by Fermi's golden 
rule, and the finite lifetimes of the energy levels are intro- 
duced phenomenologically into the density matrix equa- 



tions as relaxations. In other words, the spontaneous, 
electron tunneling induced, and stimulated emissions can 
be uniformly described by this newly proposed approach. 
The equation of motion of the density matrix is: 

dp 



-[H,p]+ TtrP + Tphp + Rp 



(1) 



Here p is the Hermitian density operator, its diagonal 
matrix element p„„ gives the probability of the molecule 
in the level n. Within the electric dipole approximation, 
the total Hamiltonian of the system H is the sum of the 
free system Hamiltonian Hq and the interaction Hamil- 
tonian of the system and the plasmon, V{r,t) = fj, • e. 
Considering the surface plasmon as a linearly polarized 
electromagnetic field, we have e{t) = eE[)Cos{ujt) . Here 
e and u are the polarization vector and the frequency of 
the surface plasmon, respectively, /x is the electric dipole 
moment operator of the molecule. Ttr is the population 
decay operator with[l5|. 



TtrP = - ^ ■ 



ajittijp - 2aijpaji + pa^iaij] (2) 



Here Ti^^j is the decay rate from the level j to the level 
2. Obj^j is a matrix with element equals to 1 and all 
other elements equal to 0. Tph is the dephasing operator 
which is introduced into the density matrix to describe 
the broadening of the emission spectrum caused by the 
surroundings. 

The electron tunneling process can only affect the pop- 
ulation of the molecule through the diagonal terms of the 
density matrix (tunnel into and out of a particular level) . 
The last term in Eq. ^ has the form, 

Rp = - y^[Rj^iajjpajj - R^^jajipaij] (3) 



Here Ri^j (Rj^i) is the electron tunneling rate from 
level i to j [j to i). Under the approximation of weak 
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coupling between the molecule and the STM tip and 
substrate, these rates can be calculated according to 
Fermi's golden rule as widely used in the rate equation 
method [l2,[ll [III . If we assume the STM tip as the left 
lead and the substrate as the right lead to the molecule. 
The tunneling rates between two relevant levels i and j 



(4) 



{Ei < Ej) could be obtained [ISH^ 



R 



with SE^ = Ej-Ei~eaV and SE'^ = Ej-Ei-e{a-l)V. 
Fij is the transition probability between the levels i and 
J. 7^ and 7^ arc the bare tunneling rates between a 
particular level of the molecule and the left and right 
leads (for the simplicity, we assume different levels have 
the same bare tunneling rates), a is the bias coupling over 
the leads, / is the Fermi distribution at the temperature 
T, fiE) = l/[l + eME/kBT)]. 

Relative to the level i in the ground state, the proba- 
bility of radiative population or depopulation of the level 
j in the excited state induced by the plasmon is depicted 
by the work of the field Wj^{t), 



(5) 



Determined by the instantaneous Rabi frequency, Wj^ {t) 
shows a sign-changing property with its positive 
{Wj'/^{t)) and negative {Wjf'{t)) parts representing re- 
spectively the absorption and the stimulated emission in- 
duced by the plasmon. Correspondingly, the work done 
by the electron tunneling process at time t can be calcu- 
lated as, 

W^M = + W^f{t) = R,^,pu - R,^rPj, (6) 

The contributions of spontaneous emissions to the to- 
tal emission spectra are described by the radiative decay 
from level j in the excited state to the level i in the 
ground state, Wjf{t) ^ —Ti^jPjj. All of these quan- 
tities can be calculated by solving the density matrix 
equation of the system, and the emission cross section 
for the transition between the level j and the level i can 
then be obtained from, 



(w^f W + w^rW + ^f (7) 



A at the right-hand side of the equation is a constant 
prefactor. The photon emission process can be described 
by a Lorcntzian function, the power spectrum is then 
obtained as, 



cr(cj) 



(8) 



Here r is the half width at half maximum (HWHM) of 
the emission line, lOji is the resonant frequency between 



the levels j and i. It is noted that the system can reach 
its steady state in several nanoseconds for all the cases 
studied here, which is several orders of magnitude shorter 
than the time of the experimental measurements. Thus 
the contribution from the transient states can be aver- 
aged out. 

We have solved the related density matrix equations for 
a TPP molecule (optimized structure shown in FiglTJa)) 
in a STM as done in the latest experiment of Dong et 
ffl/.[lj|. The three processes involved in the electrolumi- 
nescence from the molecule are schematically shown in 
FiglH The electron tunneling induced emission process 
is illustrated in Fig[T] (b) , in which electrons can tunnel 
into the LUMO and tunnel out from the HOMO when 
both orbitals are in the bias window and leave an elec- 
tronic excited state ready to fluoresce. Due to the rel- 
ative short lifetime of the vibrational excited state, the 
system will first relax to the lowest vibrational state of 
the excited state (S\) before emitting a photon accord- 
ing to the Kasha's rule. In the plasmon-assist molecular 
emission process (FiglT](c)), a plasmon can be generated 
by the inelastic tunneling of electrons from the tip to the 
metal substrate (lET process). The plasmon behaves as a 
strong coherent electromagnetic source and can excite the 
molecule resonantly into higher vibrational levels in the 
first excited state following a resonant emission, or stim- 
ulated emission, without a non-radiative relaxation 14 . 
It is noted that when the plasmonic excitation is relative 
weak, the emission will again follow the Kasha's rule. 
Spontaneous emission (fluorescence) has a similar mech- 
anism with the electron tunneling induced emission and 
it happens simultaneously with the other two emission 
processes. 

In principle, our formulation allows to include as many 
electronic states and vibrational levels as one wants to. 
However, for the relevant emission processes of the TPP 
molecule, it is enough to consider two electronic states, Sq 
and S\ (Q band), each with three vibrational levels. The 
transition probability between different vibrational levels 
is governed by the Franck-Condon (FC) factors that can 
be calculated using the standard harmonic model. The 
calculated energy gap between the two vibrational bands, 
(5(0,0) and (5(0,1), is scaled to the experimental value 
of 0.16eV and the excitation energy between two molec- 
ular states is set to the experimental value of 1.89 eV 14 1 
for the sake of a better presentation. The initial condi- 
tion for the population is po.o = 1 and pn,n = 0,n 7^ 0, 
i.e. the electrons rest at the lowest vibrational level of 
the ground state. Normally, the radiative lifetime of the 
state is about three orders of magnitude longer than those 
of non-radiative relaxations [2l|. In this work the radia- 
tive lifetime Tr is set to be 2 ns and the radiative de- 
cay rate between level j in S\ and i in 5*0 is calculated 
as ri<_j = Fji/Tr- The non-radiative decay rates be- 
tween the vibrational levels are set to be (2 ps)-^ The 
density matrix equations are solved using an iterative 
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FIG. 1: (a) The optimized structure of a tetraphenyl por- 
phyrin (TPP) molecule, (b) Electron tunneling induced emis- 
sion with the inclusion of the internal relaxation, (c) Plasmon- 
assist process, including the plasmon generation, the excita- 
tion, the stimulated emission and the spontaneous emission. 
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FIG. 2: (al) and (a2) The measured EL spectra of H2TBPP 
molecules at different bias voltages. [Adapted with permission 
from Dong et al. (Ref. 0) copyright (2004) by the American 
Physical Society.] (bl) and (b2) The calculated EL spectra of 
TPP molecules at different bias voltages without considering 
the plasmon excitation. The temperature is set to be 300 K. 



predictor-corrector method. The dephasing factor enters 
the calculated spectrum through the HWHM (r = 0.05 
eV for all the calculations in this work) of the convoluted 
Lorentzian profile (Eq. ([5])). All the calculations are per- 
formed with a bias coupling constant (a) of 0.50 and 
nearly symmetric bare tunneling rates (7^=16. 40ueV, 
^^=4.80iicV) as those used by Seldenthuis et aL[12| for 
a similar system. 

We first simulate the electroluminescence (EL) spec- 
tra of a TPP molecule in a STM induced by the elec- 
tron tunneling process only. The bias voltage depen- 
dence of EL spectra for different porphyrin molecules 
were experimentally observed by Qiu et al.^ and Dong 
et It was found that the spectral peaks under dif- 

ferent bias voltages remain unchanged with enhanced in- 
tensities at higher bias voltages. Dong et al. also re- 
ported bipolar EL spectra at ± 2.5V. In a very recent 
theoretical study by Seldenthuis et al. [l^l , the bias volt- 
age dependence of EL spectra measured by Qiu et al. 
as well as the bipolar molecular fluorescence observed 
by Dong et al. have been successfully analyzed using a 
rate equation method. We focus our attention to the 
results for the meso-tetrakis (^,5-di-terbutylphenyl) por- 
phyrin (II2TBPP) molecules [7[. Our simulated spectra 
under different bias voltages from the density matrix ap- 
proach are illustrated in Fig. [IJbl), which have largely 
reproduced the corresponding experiments as recaptured 
in Fig[2jal). The peak positions of the calculated EL 
spectra all remain the same simply due to the Kasha's 
rule. A clear enhancement of the spectral intensities at 
higher bias voltages is revealed, which is arisen from the 



increased number of vibrionic levels entering in the bias 
window. One interesting experimental observation is that 
the EL can still be observed when the excitation voltage 
is less than the energy of the emitted photon around 1.9 
eV. Such "energy forbidden" transitions were suggested 
to be resulted from either the thermally assisted tunnel- 
ing injection of holes or the shift of energy levels due to 
the charging of the molecules 0. Our simulations have 
shown that the thermally assisted electron tunneling is 
the reason because of the broad Fermi-Dirac distribution 
at BOOK. Moreover, the experimentally observed bipolar 
spectra at ± 2.50 V in Fig. [2l^a2) are also well repro- 
duced by our calculations (Fig. [2ljb2)) as a result of the 
weak couplings between the leads and the molecules. The 
intensity difference in the bipolar spectra is due to the 
slight asymmetric coupling between the molecule and the 
two leads [l^. 

It can be seen that the electron tunneling induced 
electroluminescence strictly follows the Kasha's rule, i.e. 
only the emissions from the lowest vibrational level of the 
excited state are possible. In the latest experiment, the 
EL from higher vibrational levels were observed, which 
suggests that a new mechanism must be involved [14 1. 
One of the possible candidates is certainly the plasmonic 
excitation. In order to be consistent with the experiment, 
we have considered five cases where the plasmon is res- 
onant with the transition of (0,2), (0,1), (0,0), (1,0) and 
(2,0), respectively, as shown in Fig. [3] (left). In our sim- 
ulations, these transitions correspond to the excitation 
energy of 1.57, 1.73, 1.89, 2.05 and 2.21 eV, respectively. 
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It is found from the simulations that the effect of the 
plasmonic excitation becomes visible for all the five cases 
when the field strength of the plasmon is larger than 
lO^V/m. 

The simulated EL spectra with the inclusion of both 
electron tunneling and plasmonic excitation are given in 
Fig. [3] to directly compare with its experimental coun- 
terpart. One can immediately notice that the spectral 
feature is drastically changed by the inclusion of the plas- 
monic resonant excitations. In other words, all the non- 
FC features in Fig. [31 such as new peaks, (1,0) in (h) and 
(2,0) in (j), the enhanced peaks, (0,2) in (b), (0,1) in (d) 
and (0,0) in (f), are contributed from the plasmonic exci- 
tation. We have also plotted out the contributions from 
three individual emission processes in Fig. |3] (right). It 
can be found that the electron tunneling induced emis- 
sion and the spontaneous fluorescence always follow the 
Kasha's rule and give the conventional double-peak Q 
profile. The plasmon-assist molecular emission is a stim- 
ulated process directly from the resonantly excited vi- 
brational levels. The position of such a stimulated emis- 
sion is controlled by the plasmonic resonant energy. It is 
noticed that the agreement between the theory and the 
experiment could be largely improved once the plasmon 
profile is fully considered in the simulations. 

In summary, we have presented a generalized density 
matrix formalism to successfully describe the electrolumi- 
nescence from molecules near a metal surface in a STM. 
Both the electron tunneling and local surface plasmon 
induced excitation and emission are treated on an equal 
footing. Model calculations for porphyrin derivatives 
have reproduced corresponding experimental spectra and 
nicely explained the observed unusual large variation of 
emission spectral profiles in recent experiments. 
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FIG. 3: EL spectra modified by the resonant plasmonic ex- 
citation. (Left) Measured EL spectra of TPP molecules to- 
gether with the corresponding NOP modes (shown as black 
dashed line). Adapted by permission from Macmillan Pub- 
lishers Ltd: [Nature Photonics] (Ref[i3), copyright (2009). 
(Right) Calculated EL spectra together with the plasmon- 
assist molecular emission (dashed lines), electron tunneling 
induced emission (dotted lines) and the spontaneous emis- 
sion (dash-dotted lines). The plasmon resonant energy is set 
to 1.57, L73, 1.89, 2.05 and 2.21eV, respectively, with the field 
strength of 1.0 x 10®V/m for the cases (b), (d) and (f), and 
5.0 X lO'^V/m for (h) and (j). Ah the spectra are calculated 
at +2.5 V and 80 K. 
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